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Abstract: In this letter we discuss a certain class of Hidden Valley models where the 
dynamics in the hidden sector are close to a strongly coupled conformal fixed point. We 
show that these models can explain the excess in cosmic ray muon events with high muon 
multiplicities that has been reported from cosmic ray studies with the ALEPH and DELPHI 
detectors. We further point out that these models can also at least partially account for the 
excess in multi-muon events that was found by the CDF experiment. Finally we discuss 
possible signatures of these models at the LHC. 
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1. Introduction 

Recently the CDF collaboration presented a study of multi-muon events in the bb search 
channel [1]. They considered events with at least two central muons (|?7| < 0.7), a transverse 
momentum of pT ^ 3 GeV, and an invariant mass of the muon pair greater than 5 GeV. 
Events where both triggered muons are produced within the beam pipe are well understood 
and the measured event rate is in agreement with the theoretical prediction. 

Once the requirement that the two muons are produced in the beam pipe is removed, 
the observed number of events lies above the theoretical prediction - about 150000 events 
remain unexplained and have been termed the ghost sample in [1]. They further observed 
that within the ghost sample, the number of additional muons within a cos 9 > 0.8 cone 
around a primary muon is unexpectedly large, more than four times what one would expect 
for an ordinary standard model event. 

CDF is not the only experiment where an excess of multi muon final states in hadronic 
collisions was found and may hint to new physics. In [2] the DELPHI collaboration at LEP 
studied multi-muon events coming from cosmic ray showers. They compared the number 
of events with a certain muon multiplicity with the prediction of a Standard Model Monte 
Carlo simulation. They found that while the simulation agrees with the number of events 
for low and medium muon multiplicities, the number of events with the highest muon mul- 
tiplicities, (A'^ > 45), clearly shows an excess over the theoretical prediction. Depending 
on the model it is an excess of up to 3<t compared to the Standard Model. 
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The findings of tlic DELPHI study [2] nicely confirm the result of an earlier study by the 
ALEPH collaboration, which observed [3] an increase of high multiplicity muon events com- 
pared to the Monte Carlo prediction. Although they can explain the integral multiplicity 
distribution for multiplicities of < 60 they fail to explain the highest multiplicities. 

Hidden valley models [4] are a class of hidden sector models with strong dynamics and 
a spectrum light enough to be observable at current or future collider experiments. Their 
phenomenology, in particular at hadron colliders, has been studied in [5-10]. In [11] it has 
been pointed out that the CDF multi-muon excess could be explained by a Hidden Valley 
model. 

In this paper we analyse a certain Hidden Valley scenario with emphasis on its phe- 
nomenology in cosmic ray showers. We show that the excess of multi-muon events can be 
explained by the production of a heavy particle that decays into the hidden sector. We 
further argue that the assumptions needed to account for the cosmic muon excess also offer 
the possibility that this model explains the multi-muon excess reported by CDF. 

This paper is organised as follows. After discussing the muon excess in cosmic ray 
showers in the next section, we show how conformal Hidden Valley models can explain 
that excess and argue that it could also explain the CDF anomaly. In section 4 we discuss 
possible LHC signals of these models before we conclude in section 5. 

2. Cosmic Rays and Muons 

Cosmic rays originate from extraterrestrial sources that seem to be uniformly distributed. 
These particles can be electromagnetic particles such as electrons and photons as well as 
strongly interacting particles such as protons and nuclei. For energies above ~ 1 GeV the 
electromagnetic component is negligible so that the usual assumption is that the incident 
particle is a proton or a nucleus. The most stable nucleus and the heaviest that can be 
produced by nuclear fusion is iron. Therefore one assumes that the cosmic rays consist of 
a mixture of nuclei from single protons to iron. A more detailed discussion of cosmic ray 
showers and simulations of the muon density profile can be found in [12]. 

Detected muons from cosmic ray showers are mainly decay products from hadrons pro- 
duced after the initial collision of a cosmic particle with a nucleus of the atmosphere. The 
setup of the DELPHI experiment allows to observe events with cosmic ray energies between 
10^^ — 10^^ eV. These high energies assure that the muon bundle is strongly boosted and 
this restricts the radius of the cone of the decay products. The hadronic showering and 
the subsequent decay into muons can be simulated using Monte Carlo methods. The DEL- 
PHI collaboration used QGSJETOl [13] to simulate the initial collision and the COSIKA 
program [14] to model the atmospheric shower. They have varied the center of the shower 
within a certain radius around the detector. Their result was that the optimal radius 
is around R = 200 m. If the radius is smaller, one looses events with low multiplicities 
whereas if the radius is larger one needs larger data samples to also get events with high 
multiplicities. For that value of the radius one looses less than 0.5% of the events with the 
lowest multiplicity of = 4. This leads to the conclusion that the muons in a cosmic 
shower are typically concentrated in a radius of 200 m from the center. 
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From these simulations one can nicely see that for very low multiplicities the data fits 
the model with a proton-induced shower whereas for multiplicities with 5 < A^^ < 50 the 
data is well described by an iron induced shower. However for multiplicities with A?^ > 50 
the simulation fails to describe the data. 

A similar characteristic behavior has been found by analysing cosmic ray data of the 
ALEPH experiement [3], where up to a multiplicity of A^^ < 50 data and simulation agree 
well. But as for the DELPHI detector there is an excess of muons for higher multiplicities 
that cannot be described by the simulation. 

3. Hidden Valleys and Muons 
3.1 Basics 

Hidden sectors have been known in particle physics for quite some time. The basic idea is 
that besides the standard model sector there exists another hidden or dark sector in the 
theory that only communicates with the SM via some messenger field X. Usually X is 
assumed to be quite heavy or couples only very weakly to the SM, and the non-observation 
of the hidden sector is attributed to this fact. This type of construction is most popular in 
models of supersymmetry breaking, where supersymmetry is broken in the hidden sector 
and communicated to the visible sector via the messenger fields. In most of these models 
the hidden sector is assumed to be heavy and its phenomenology is neglected. 

Only recently the possibility that the hidden sector could be relatively light has been 
considered. In the so-called Hidden Valley models [4] a new confining gauge group Gy 
is added to the Standard Model that communicates with the SM only via a heavy Z'. 
New fermions, u-quarks, that are charged only under are confined and form i;-hadrons. 
The properties of the v-hadrons are highly model-dependent, but the authors of [4] have 
identified some typical properties. These include long lived i;-hadrons that decay back to 
SM particles with a displaced vertex and possibly also stable v-hadrons that provide dark 
matter candidates. 

Further support for Hidden Valley models with non-negligible phenomenology comes 
from the recent experimental results of the PAMELA experiment that reports an excess 
of positrons in cosmic radiation that could originate from dark matter annihilation in the 
galactic halo [16]. The authors of [17] have shown that this can nicely be explained in 
terms of a hidden sector model of dark matter. In their model, dark matter is part of a 
dark sector with a gauge group Gdark- The corresponding dark gauge boson mixes weakly 
into the photon and by this mixing the sector couples to the standard model. The weak 
mixing allows the dark gauge boson as light as 1 GeV or even less, which is required to 
successfully explain the positron excess in terms of dark matter annihilation. 

The simple Hidden Valley model [4] described above assumes QCD-like dynamics. 
Hadronization will therefore proceed similar to hadronization in QCD, so after a collision 
in most cases two v-hadron jets will form, and if the lifetimes of the v-hadrons are not 
too large, will yield events similar to those coming from a qq ^ qq scattering. However, 
the strongly coupled sector in the hidden valley could behave completely different from 
what is known from QCD. This will for example be the case if the Gy gauge theory is not 
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asymptotically free but confining in the UV instead, or, if the theory lives close to a strongly 
coupled fixed point. These cases have been studied in [18], and it was argued that the effects 
could be quite spectacular. The reason is that at large coupling (or more precicely, at large 
t'Hooft coupling g^N , where is the number of colors of the hidden gauge group G^) 
the parton shower becomes much more efficient and thus all quarks and gluons rapidly 
become soft. Instead of two or three hard jets this then yields an approximately spherical 
distribution of soft u-hadrons that decay back to SM particles with an almost uniform 
angular distribution. 

3.2 Hidden sector model 

For our purposes we will assume a hidden sector that has the properties described in the 
previous paragraph. We assume that a heavy resonance X with a mass mx ^ 200 GeV 
can be produced in pp collisions, either alone or accompained by a quark or a quark pair, 
with a quite large cross section. The resonance has a large (> 99%) branching fraction 
into w-quark pairs that immediately shower and yield a nearly spherical distribution of v- 
hadrons. We further assume that the main decay mode of v-hadrons back to SM particles 
is into leptons, in particular muons. This could be realized for example by a mechanism 
similar to the one presented in [17]. For definiteness we assume that the number of SM 
particles from one X decay of the order of 100, and that about half of these particles are 
muons. This number will of course fluctuate largely from event to event. 

One immediate question is of course whether these events should have been seen in 
other experiments, in particular at colliders. The mass rax is chosen such that a production 
of this particle at LEP is highly unlikely, if not impossible. And it is not clear whether 
these events could have been detected by the LEP experiments even if mx would be much 
lower. The reason is that only a few if any of the produced SM particles will have energies 
above one GeV, and electromagnetically charged particles may not have been sufficiently 
energetic to be detected. In this context, it is interesting to note that DELPHI also observes 
a soft-photon excess in hadronic events. The same argument in principle applies to other 



experiments, however in section 3.4 we argue that CDF might be able to see some of these 



muons and that they even could account for a part of the CDF muon anomaly. 
3.3 Explaining the cosmic ray muon excess 

At collision energies above mx the X resonance is copiously produced and decays into N 
muons and some additional SM particles. The number may vary from event to event, but 
is assumed to be of order ~ 50. The actual collision energy between the partons inside 
the proton is on average much lower than the center of mass energy of the proton-proton 
collision itself. To obtain high enough rates, we thus require Ecm ^ lOmx- For cosmic 
rays consisting of protons this translates into a lower bound on the energy of approximately 
Ep > 10^ GeV, which coincides with the energy range that accounts for most of the events 
with > 45. For cosmic rays consisting of iron nuclei the bound is much higher, thus 
we require that at least a large fraction of cosmic rays consists of single protons or light 
nuclei. The chemical composition of cosmic rays has been studied by the MACRO and 
EAS-TOP collaboration [15]. They found that cosmic rays can not consist of a single 
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Figure 1: Maximal distance Rmax of the muon point of impact from the shower centre, as a 
function of the muon energy in the rest frame of X, shown for three different boost factors and a 
distance L = 10 km from the surface. 

particle species and that a mixture of light and heavy nuclei provides the best fit to the 
data, which supports our assumption. 

In its center of mass frame the X resonance decays almost isotropically. The energy is 
distributed more or less equally between the 0(50) muons, thus the energy per muon lies 
in the range ~ 1 — 5 GeV. To be able to explain the excess in events with a large muon 
multiplicity, we must require that most of these muons hit the DELPHI detector, i.e. after 
travelling 0{W'^m) through the atmosphere they should not deviate more than 5 meters 
from the track of the incoming cosmic particle. 

To be more definite, we consider the following setup: The cosmic ray propagates in 
the z-direction with 4-vector pi^^ = {El,0, 0, El), neglecting its mass, and hits a proton at 
rest with mass rup. We denote quantities in the laboratory (earth) frame with a subscript 
L, and quantities with respect to the center of mass of the initial collision with a subscript 
CM. The boost factor into the center of mass frame of the collision is approximately given 



This boost factor is of order 10^ — 10"^ for the relevant energy range. For simplicity we first 
assume that the X resonance is produced at rest. The distance R between the point where 
the muon hits the ground and the centre of the cosmic ray shower then only depends on 
the energy of the muon and on the boost factor, i.e. on the energy of the incoming cosmic 
ray. 

In this scenario the boost factors are large enough so that all muons hit the surface 
within a distance Rmax from the shower centre. The dependence of Rmax on the muon 
energy for different boost factors is shown in figure |l|. 



by 
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Figure 2: The first figure shows the radial muon flux density for a muon rest frame energy 
= 2 GcV for three different boost factors and L = 10 km, normalised to j{0) = 1. The second 
figure shows the integrated distribution, i.e. the fraction of all muons that hit within a radius R of 
the shower centre, also for L = 10 km. 



To determine the fraction of muons that actually hit the detector we need to know the 
radial flux density distribution j{R) for a given boost factor and muon energy. Figure ^a 
shows this quantity for a muon energy of = 2 GeV and different boost factors. The 
distributions are almost identical for different muon energies. The integrated distribution, 
i.e. the fraction of muons that hit within a radius R of the centre is shown in figure Qb. The 
calculation of the radial and integrated distributions is detailed in the appendix, where we 
also give the explicit formulas. 

From figure |2| we can easily see that for a boost factor 7 = 3000 more than two thirds 
of all muons hit within a radius of 5 meters from the centre, while for a boost factor larger 
than 7000 more than 90% of all muons are within that range. A boost factor of 7 = 3000 
corresponds to an energy slightly above 10^ GeV. 

The analysis of [2] has shown that that the majority of the events with Nfj_ > 45 comes 
from cosmic rays with an energy of 10^ — 10^ GeV and with the shower centre at or very 
close to the DELPHI detector. These are precisely the events where at least a large fraction 
of the muons coming from the X decay would also hit the detector, since we have shown 
above that for this energy range most of the additional muons come within a radius of 5 
meters from the centre. 

We still need to estimate the fraction of events that contain an additional X particle 
such that the excess in multi-muon events can be explained. In the DELPHI event sample 
the number of events yielding 45 or more muons is around 350, the number of events with 
Xfj, > 30 is 1065 from table 1 in [2]. We thus assume that at least 500 events fall into the 
category described above, i.e. events with an initial energy > W''' GeV that hit the 
detector almost centrally. The number of observed events with Nfj_ > 70 and Nfj_ > 100 can 
be taken from the same table, while we take the number of expected events for the same 
region from figure 14a in [2]. The actual numbers are shown in Tab. ||. Assume that a 
fraction k of the above events contains an X resonance, and that in these cases on average 
40 additional muons hit the detector. One can then see that with a k of order 5% — 10% 
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experiment 


proton, theory 


iron, theory 


> 70 


78 


26 


56 


> 100 


24 


6-7 


13 



Table 1: Observed and expected number of events at DELPHI [2] with more than 70 muons and 
more than 100 muons respectively. 

the X resonance decays can account for the additional high muon-multiplicity events. 

Of course this is only a rough estimate. From the energy dependence of high energy 
hadron collisions, we know that k will not be a constant but most likely will increase for 
increasing cosmic ray energies. This could further enhance the number of events in the 
highest multiplicity range, > 100, which seems to be indicated by the data. In the end 
a full simulation of cosmic ray showers with a production of an X resonance is needed to 
see whether this model agrees completely with the data. 

3.4 The CDF multi-muon excess 

In Ref. [1] the CDF collaboration studied di-muon and multi-muon events in the bb search 
channel. The sample consists of events with at least two muons with an invariant mass 
5 GeV < rrinn < 80 GeV and transverse momentum pt > 3 GcV. They further split the 
sample into two parts, one part where both muons are produced within the beam pipe, 
and the other part where at least one muon originates from outside the beam pipe. The 
multiplicity of events where both muons are produced inside the beam pipe are perfectly 
described by QCD production [1]. Once the full sample is studied, approximately 150000 
events remain unexplained. These events are called ghost sample in [1] . About half of the 
ghost sample can be accounted for by fake muon signals and other production mechanisms. 

Additional muons have then been searched in a cone of an angle 9 < 36.8°, {cos 9 > 0.8) 
around any of the initial muons in the ghost sample. It was found that the average number 
of additional muons in the ghost sample are unexpectedly large, at least four times of what 
one would expect for an ordinary QCD event. 

The question now is if these events can be explained by a hidden sector model along the 
lines of section 3.2. We assume that the X resonance is also produced in the pp collision at 
Tcvatron and decays into 0(50) muons. As the hidden sector particle is heavy, of the order 
of 100 — 200 GeV, there is not a large boost factor between the ccntrc-of-mass frame and 
the laboratory frame. Therefore the particle will decay rather isotropically. It is further 
reasonable to assume that most of these muons will have a somewhat displaced vertex since 
they originate from a showering process in the hidden sector. 

It is likely that one pair of muons from an X decay has an invariant mass above the 
required 5 GeV. This determines the number of X particles that must be produced in 
order to explain the ghost sample to be around 70000. This number is quite large but not 
completely unreasonable, as has been shown for example in [19]. Due to its exotic decay 
mode, we do not expect this process to be excluded by any existing new physics search. 
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The next question then is the multiphcity of additional muons in these events. A cone 
with the angle of cos 9 > 0.8 covers 10% of the solid angle. Summing both initial muons 
then cover 20%. We therefore expect 0(10) additional muons in the cones for such a decay. 
However most of these muons will not pass the required pt ^ 3 GeV, since, as discussed 
above, the average energy per muon is assumed to be around 2 GeV, but could also be 
much less. Since we do not know the energy distribution of the muons we can only give 
a rough estimate here. Assume that 10% of the muons passes the pT cut. This gives us 
0(5) muons. Two of these will be the initial muons that make up the ghost sample. Of 
the remaining muons, only about 20% will be in a cone around the initial muons. It is 
therefore not unlikely that the probability for a third muon is between 10% and 30%, and 
further the probability for a fourth muon is again ~ 20% of the three muon probability. 
Thus at least qualitatively this signal is similar to the one found by the CDF analysis. 

We do not dare to draw a final conclusion here, since too many unknowns remain. 
One could for example consider processes where the resonance is produced along with a 
recoiling particle and a sizeable boost factor. This could largely increase the number of 
events with very high muon multiplicities, for > 200 GeV these events would however 
suffer from strong kinematical suppression at Tevatron energies. We will present a more 
detailed study of such events at hadron colliders elsewhere. 

4. Multi-Muon events at LHC 

For the Tevatron the centre of mass energy of 1.96 TeV justifies the assumption that the 
hidden sector particle is produced nearly at rest and therefore rest frame and laboratory 
frame are practically identical. This statement however no longer holds for LHC energies. 
Here it is much more likely that the X resonance is produced with a moderate to large 
boost relative to the detector. Especially, the available collision energy should be sufficient 
to produce the X resonance recoiling against an ordinary QCD jet. This will lead to events 
as depicted in figure ^, where a large number of muons with energies around 10 GeV are 
radiated in the same direction, forming a muon jet [17]. 

The angular distribution of muons depends mainly on the boost factor. This also 
determines the average number of muons that will hit in a cone around the boost direction. 
In figure ^ we show the fully differential angular distribution of the muons for different boost 
factors. The angle 9 always denotes the polar angle with respect to the boost direction 
and with = lying in direction of the boost. For large boosts one observes that most 
muons will end up in a rather narrow cone. The muon energies are around 5 — 10 GeV 
for smaller boost factors and up to 50 GeV for the most energetic muons and 7 = 4. The 
signal therefore is a large muon jet with single particle energies between 5 GeV and 50 GeV 
and a small number of softer muons into other directions that might remain undetected 
due to their lower energy. 

Events where the multi-muon final state is boosted to sufficiently high transverse mo- 
menta to be fully contained in the detector could moreover offer the very exciting prospect 
of determining the dynamics of the u-parton cascade, thereby obtaining information on the 
hidden gauge group. To do so, one would reconstruct the momenta of all produced muons, 
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Figure 3: For LHC energies the hidden sector particle can be boosted in the laboratory frame and 
therefore also leave a characteristic signature in the detector depending on the boost factor (see 
figure 3). 




Figure 4: Angular distributions of muons in the laboratory frame, for different boost factors, as a 
function of the angle 9 with respect to the boost direction, p is the momentum of the recoiling jet. 
An unboosted decay would produce a straight line. Distributions arc normalised to unity. 



boost them back into the rest frame of X, and measure the classical event shape variables 
(thrust, C-parameter, oblateness, broadening) from these rest frame momenta. 

5. Conclusion 

Hidden Valley models with conformal dynamics offer a vast and spectacular phenomenol- 
ogy, that is largely unexplored up to now. In this letter we have introduced a variant of 
such a model and have shown that it can explain the excess of high multiplicity muon 
events in cosmic ray showers. We have further shown that the multi-muon events recently 
observed by the CDF collaboration could in principle originate from the decay of such a 
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hidden sector resonance. A more rigorous analysis is required to find out whether both 
these signals indeed have a common and whether a conformal hidden sector will turn out to 
be a viable explanation. Our work is the first step into this direction. Other explanations 
for multi-muon production in cosmic ray events have been given in terms of strangelets [21] 
and in terms of anomalous multi VF-boson production [22] . We have seen here that hidden 
sector models typically predict the production of soft muons in the CM frame of the cosmic 
ray interaction, and that they predict a boost-dependent characteristic radial profile of the 
muon flux. Future measurements of muon production in cosmic rays with detectors at 
the surface or at shallow depth could provide further insight into the dynamics of these 
anomalous multi-muon events. 

Alternative models have been proposed to explain the CDF anomaly [16,20], based on 
cascade decays of new particles into 8 or 10 tau leptons. While these models could turn out 
to be more successful in explaining the CDF signal, they can only marginally contribute 
to the muon multiplicity in cosmic ray events and therefore cannot explain the excess seen 
at DELPHI and ALEPH. 

Finally we pointed out the possibly exciting signatures that such a model will produce 
at the LHC. 
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A. Boosted Muon Densities 

Assume that the muon is emitted with an energy of in the CM frame of the decaying 
X. Then its momentum is given by 



p'^^ = {E^,psme,0,pcos6) , p=^El-m'^, (A.l) 

where 6 is the angle in the CM frame under which the muon is emitted. In the laboratory 
system the momentum is given by 

p"^ = {^{E^-Pp cose), p sine A 7{p cose- pE^)f, p=^El- w?. (A.2) 

where 7 is the boost factor from the CM frame into the laboratory frame, and (3 = 
^yi — 7~^. Typical boost factors for high energy cosmic ray events range from 7 = 10^ to 

7 = lOl 

Assume the event takes place at a height h from the surface. For large enough boost 
factors all muons will hit the surface. The distance of the point of impact from the shower 
center depends, for fixed energy, only on the angle e under which the muon was emitted, 
and is given by 

m = h^M = h , (A.3) 
v'iiQ) 7(/3f-cos0) ^ ^ 
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For small boost factors or for massless particles this expression has a singularity at cos = 
(3E/p that determines an angle Omin- Muons with 9 < 9min don't reach the surface but 
are radiated back into the universe. For massive particles and large enough boost factors, 



p 

that is given by 



such that {(3— — cosO) > for all 6, the expression has a maximum Rmax at an angle Or, 



Omax = arccos {^-^ ' (^-4) 

and Rmax — Rifimax^ 

We now calculate the flux density as a function of the distance R. It is given by 
dN 

j{R) = — with dN = 2TTsm9de, dA = 2TrRdR. (A.5) 



Inserting dN and dA we obtain 



For 7 = we have R = h tan 9 and a straightforward calculation yields 

j{R) = ^cos^e{R). (A.7) 

For 7 7^ the function dO/dR can be determined from equation ( |A.3| ). The final 
expression is 

sin 9 7 / cos 9 sin^ 9 \ 

^^""^ = ^ ^ X I (/3f - cos^) " (/3f - cos 0)2 j ' ^^-'^ 



where 9 = 9{R) has to be determined from the inverse of eqn. ( A.3| ). The function R{9) is 



not monotonic over the full range of 9. One has to distinct the two cases 9 < 9 max and 
G > Omax and add both contributions to obtain the full density profile. As it turns out for 
large boost factors the contribution from < 9max can be neglected. 
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